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Fluctuations around a Bose-Einstein condensate can be described by means of Bogolubov the- 
ory leading to the notion of quasiparticle and antiquasiparticle familiar to non-relativistic condensed 
matter practitioners. On the other hand, we already know that these perturbations evolve according 
to a relativistic Klein-Gordon equation in the long wavelength approximation. For shorter wave- 
lengths, we show that this equation acquires nontrivial corrections which modify the Klein-Gordon 
product. In this approach, quasiparticles can also be defined (up to the standard ambiguities due 
to observer-dependence). We demonstrate that — in the low as well as in the high energy regimes — 
both concepts of quasiparticle are actually the same, regardless of the formalism (Bogolubov or 
Klein-Gordon) used to describe them. These results also apply to any barotropic, inviscid, irrota- 
tional fluid, with or without quantum potential. Finally, we illustrate how the quantization of these 
systems of quasiparticles proceeds by analyzing a stationary configuration containing an acoustic 
horizon. We show that there are several possible choices of a regular vacuum state, including a 
regular generalization of the Boulware vacuum. Issues such us Hawking radiation crucially depend 
on this vacuum choice. 



PACS numbers: 04.62. -fv, 03.75.Kk, 04.70.Dy 

I. INTRODUCTION 

It has for quite some time been understood that 
the propagation of sound waves in inviscid irrotational 
barotropic fluids can, under rather general circum- 
stances, be effectively described by relativistic curved- 
spacetime metrics, see e.g. [T]. This led to the obser- 
vation that it should, at least in principle, be possible 
to achieve sonic or acoustic black hole configurations, 
and thus study certain types of high-energy effects on 
black hole physics by analogy This gravitational 

analogy was established on firm foot about a decade ago 
in |3] . Since then, the field of analogue gravity has grown 
into a mature and well-established research programme 
within the gravity community [4j. It offers both the ex- 
citing prospect of bringing black hole physics, in partic- 
ular Hawking radiation, within the reach of experimen- 
tation, as well as a refreshing conceptual take on some 
long-standing problems of quantum gravity, such as its 
relation to dark energy [5j or the avoidance of black hole 
singularities in the gravitational collapse of ultra-heavy 
bodies [fj!. 

The quest for a useful background fluid in which to con- 
duct experiments of analogue gravity leads in the direc- 
tion of superfluids. Indeed, superfluids have a vanishing 
viscosity, and some superfluids can be made extremely 
pure. Among the possible candidates, Bose-Einstein con- 
densates [3 [5] have the comparative advantage of be- 
ing conceptually well understood, and relatively simple 
to describe theoretically and manipulate experimentally. 
Acoustic black holes have been recently reported for the 



first time in Bose-Einstein condensates [9], and there is 
good hope that an experimental detection of analogue 
Hawking radiation, one of the main current promises of 
the analogue gravity programme, will be achievable in 
the near future [10] . 

Another motivation for studying Bose-Einstein con- 
densates in the context of analogue gravity is that they 
provide a real system in which high-energy modifications 
to the relativistic dispersion relations arise. This shows 
that there exist concrete examples in nature in which 
Lorentz invariance is realized as a low-energy effective 
symmetry, broken at high energy. Such a scenario could 
be the case for the local Lorentz invariance of general 
relativity as well. In this manner, Bose-Einstein con- 
densates might offer an interesting model for quantum 
gravity phenomenology, see e.g. [TTlfT^ . 

Our first aim in this paper is to further elaborate and 
consolidate the theoretical framework for the gravita- 
tional analogy in Bose-Einstein condensates by taking 
a look at the different inner products and related cre- 
ation and annihilation variables that can be introduced 
depending on the point of view. Quantum sound in Bose- 
Einstein condensates can, on the one hand, be analyzed 
within the Bogolubov formalism by directly perturbing 
the Gross-Pitaevskii equation. On the other hand, the 
phase perturbations of the condensate obey a modified 
Klein-Gordon equation and a corresponding quantiza- 
tion can be carried out. Remarkably, both procedures 
give rise to the same quantum theory. This allows us 
to establish a deep conceptual connection between both 
formalisms, the first one being inherently non-relativistic 
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while the second is relativistic, up to corrections which 
are vanishingly small for long wavelengths. Both proce- 
dures are known to be equivalent in this long-wavelength 
acoustic limit, see e.g. [13 . Here we discuss in detail how 
the usual Klein-Gordon field is distorted after the acous- 
tic approximation is broken. The step-by-step analysis 
that we perform shows that the equivalence between the 
Bogolubov approach and this generalized Klein-Gordon 
formalism persists well beyond the limit of validity of the 
acoustic approximation. 

Once an inner product has been defined, one can pro- 
ceed with the quantization of the system of quasiparticles 
in the standard Fock manner. An orthonormal and com- 
plete set of positive norm modes is needed to expand the 
field operator. Finding such a complete set of positive 
norm modes amounts to defining quasiparticle creation 
and annihilation operators and a vacuum state. This con- 
struction is not unique: one can choose several distinct 
sets of modes and vacuum states which give place to dif- 
ferent quasiparticle notions. To illustrate this procedure 
we will consider a stationary one-dimensional configura- 
tion possessing an acoustic black hole horizon. The sec- 
ond aim of this paper is precisely to show that for this 
configuration there exist several regular vacuum states 
with a specific interpretation, including a regular gener- 
alization of the Boulware state for a relativistic field in a 
black hole geometry. In this way, we show that the free- 
dom in choosing a vacuum state is larger in dispersive 
theories than in relativistic theories. This vacuum choice 
has crucial importance in issues such as the presence or 
not of Hawking radiation. 

This paper is organized as follows. Section II is de- 
voted to the Bogolubov approach and the definition of 
the appropriate inner product. A mode analysis is also 
performed. Section III introduces the hydrodynamic rep- 
resentation that leads to a Klein-Gordon equation in 
the long wavelength regime and to its generalization for 
all wavelengths. The corresponding generalized Klein- 
Gordon product is also introduced. In Section IV, we 
show that both formalisms are actually equivalent and 
lead to the same concept of positive and negative norm 
solutions. After defining the appropriate inner product. 
Section V is devoted to a discussion of the different vac- 
uum state choices. We focus on a configuration of partic- 
ular relevance for analogue gravity experiments: a one- 
dimensional stationary acoustic black hole. We conclude 
with some final comments in Section VI. 



II. BOGOLUBOV APPROACH 

Let us consider a condensed dilute gas of interacting 
bosons described in terms of quantum field operators 
■f/) and that annihilate and create particles (see e.g. 
Refs. [m m]). The operator tjj can be separated into 
two parts: a macroscopic wave function i/jq describing 
the actual Bose-Einstein condensate and a quantum field 
operator describing perturbations around the conden- 



sate. 

The order parameter ipQ satisfies the Gross-Pitaevskii 
equation 

iHdt4'0 - f-;^V2 + + g\i;l\) V-o , (1) 
\ 2m / 

where T4xt is an external potential, m the atomic mass 
and g the atomic interaction constant (proportional to 
the s-wave scattering length). We will also use the 
Madelung representation 

i'o = V^e'^"/" , (2) 

in terms of the number density of atoms uq and the phase 
Oq of the condensate, which defines the flow velocity po- 
tential (such that the flow velocity is v = \/9q /m) . Then 
Eq. ([I]) translates into 

dtm + w ■{mweo)/m = ,(3) 

dtOo + :^ (V0o)' + Kxt + gno - |^ = , (4) 

2m 2m ^Juq 

which are the continuity equation and the Bernoulli equa- 
tion plus a quantum potential term, respectively. 

A. Bogolubov equation 

The quantum perturbation field (/> satisfies the Bogol- 
ubov equation (see e.g. Refs. [HI [T5] l 

ihdt'^ ^H^ + TOc^e^*^"/''^!' , (5) 

where — gno/m is the square of the local comoving 
speed of sound and H is the operator 

n = -T^V^ + -mv2 - dt9o + mc" . 6 

2m 2m c 2 

The commutation relation for the atomic creation and 
annihilation operators 

[^(x,i),^^(x',t)] -<5(x-x') (7) 

translates into the commutation relation 

[0(x,i),0^(x',i)]=5(x-x') (8) 

for the perturbation field operators, which create or an- 
nihilate atoms in the noncondensed part and correspond- 
ingly annihilate or create them in the condensed phase. 
In other words, ((> moves an atom from the noncondensed 
part to the condensate, and vice versa for 0^. This com- 
mutation relation is valid for condensed systems in which 
the number of noncondensed atoms is very small com- 
pared to the number of condensed ones. 

Note that the Bogolubov equation ([5| could also have 
been obtained by expanding to first order in (p the 
Gross-Pitaevskii equation (IT]) for the order parameter 
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_y^giSo/?i _|_ other words, the classical perturba- 
tion (f) of the mean-field wave function of the condensate 
satisfies exactly the same evolution equation as the quan- 
tum fluctuations (j) around the condensed phase. From 
now on, we will therefore drop the hat from the operator 
(j), unless necessary. 

It is important to note that, although this equation can 
be obtained by 'linearizing' the Gross-Pitaevskii equa- 
tion, it is a complex equation for a genuinely complex 
field and is therefore non-linear: If is a solution, then 
in general a(j) is not (unless a is real). Therefore, we 
cannot directly perform a mode expansion to find the 
general solution. There exists a procedure [M f lTS ] which 
allows to overcome this problem by enlarging the space 
in which we look for solutions to Eq. ^ and to define an 
inner product in this enlarged space. 

With this aim let us introduce the spinor field 



1 



subject to the evolution equation 

where A4 is the operator 
H 



(9) 



(10) 



M 



-n 



- • (11) 



In this equation, a± = (cTx ± iffy)/ 2 and <Tx,y,z are the 
Pauli matrices. 

This equation is now linear, i.e., if $ is a solution, then 
so is for any complex constant a. The solutions to 
the Bogolubov equation ([5| are obtained by restricting 
the solutions of Eq. ([To| by the condition 



I.e. 



(12) 



B. Bogolubov inner product 



Taking into account that (j) has been defined as an op- 
erator that annihilates atoms in the noncondensed part of 
the gas (and hence creates them in the condensed part), 
the expectation value of will provide the number of 
noncondensed atoms iVi (under the assumption that this 
number is small compared to the total number of atoms, 
as mentioned above). This condition translates into the 
following normalization condition for the spinor $: 



in the standard inner product 



(13) 



($!$') d^x$^$' = d^x (</)*(/.' +^*^') , (14) 



where D is the number of spatial dimensions under con- 
sideration. Since there is a continuous exchange of atoms 
between the condensed and noncondensed phases, this 
norm will not be conserved in time. Indeed, this fact is 



already encoded in the evolution equation ( 10 1. Actually, 
it is straightforward to see that 



(15) 



However, the operator A4 is not selfadjoint in the posi- 
tive definite inner product ( 14 1, but satisfies the following 
properties: 



(JxMcFx = -M* 



(16) 



and hence this inner product is not preserved in the evo- 
lution, as we had advanced. 

These are the relevant properties of A^. In fact, in view 
of the last equality, we can introduce a "Bogolubov" inner 
product 



($i$')b = y"d^x$v,$' , 



(17) 



in which M is selfadjoint. Indeed, it is straightforward 
to check that 



(18) 



Therefore we see that the price to pay for making the 
evolution operator Ai selfadjoint is the introduction of an 
inner product (•|-)b which is not positive definite. Indeed, 
this Bogolubov inner product has the following proper- 
ties: 

• It is conserved in the evolution of the lab time t. 



• It is hermitian, i.e., ($|(f>' 



• It is antilinear in the first argument and linear in 
the second, i.e., for any complex number a, 

(a$|$')B = a*($|*')B , ($|a$')B = (19) 

• It is not positive definite, since it satisfies 

(f7,$*|a,$'*)B = -($'|$)b . (20) 

Finally, note that the physical solutions, i.e. those that 
satisfy ax^* = $ because of condition (12), have zero 
norm, as can easily be seen from Eq. ( 20 1 . 



C. Mode expansion 

The evolution operator M is selfadjoint in a non- 
positive-definitc inner product and therefore it may have 
complex eigenvalues. We will assume that the condensate 
is stable, which implies that genuinely complex frequen- 
cies cannot be present. 
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In view of the properties (16), it is easy to see that, if 



1 (uk 



(21) 



\/2 

is an eigenspinor of with eigenvalue oj/c, i.e., if 

MUk = ^fciJfe , (22) 

then 

1 



72 V"fc 



(23) 



is an eigenspinor of M with eigenvahie —LOk- Besides, 
(JzUk is an eigenvector of with eigenvalue Wfc. Fur- 
thermore, the modes Uk and are orthogonal (and can 
be chosen orthonormal) in the Bogolubov inner product: 



{Uk\Ui)^ = 2] (ul^i ~ vlvi) ^ Ski , (24) 
{Uk\VnB -'11^''^ - = ' (25) 

{V:\VnB A^'xivkvt -uuu*i)^~5ki . (26) 



Any spinor $, solution to Eq. (10), can be expanded in 
this basis: 



and its norm is given by 



($|$)3 = ^(|a,p-|fe,n 



(27) 



(28) 



We therefore again see that the physical solutions — those 
satisfying the condition ( 12 ) — have zero norm, since they 
satisfy = 6^. Note that the modes themselves are not 
physical, not only because they may be generalized eigen- 
vectors normalized to the Dirac delta, but also because 
in general they do not satisfy Eq. ( 12 ): OxU^ = ^ Uk, 
as we have seen. 

Finally, the number of field degrees of freedom carried 
by the spinor <I> is just two. Indeed, two complex (four 
real) functions are needed at an initial time to obtain 
the value of $ at any other time. Condition (12), which 



ensures the physical nature of the configuration, reduces 
this number to one (two real initial functions), which 
is precisely the number of field degrees of freedom of a 
real relativistic scalar field. In the next section we will 
actually describe the condensate perturbations as a real 
scalar field satisfying (in the appropriate limit of long 
wavelengths) a relativistic wave equation. 



III. KLEIN-GORDON APPROACH 



Bernoulli ^ equations around a background condensate 
characterized by no and ^o- Let us introduce the density 
hi and phase Oi perturbations: 



n = no + g ^rii , 9 = 9o + Oi ■ 

These perturbations obey the equations 

dtfii + V • (niv + c^Veii) = , 
at6'i-^v-V6li + (l-e)ni =0 

where Q is the operator 

e-Je^v[c2.v(c-M], 



(29) 



(30) 
(31) 



(32) 



^ = h/{mc) is the healing length, and the 7k- stands for 
the argument upon which Q acts. 

Note that the perturbation fields 0, introduced in the 
previous sections are related to the complexified density 
hi and phase 9i perturbations in the following way: 



-iOo/h 



1 1^ 

^ 7^ni+t-9i 



(33) 



in terms of which 



hi = Cy/gm{e 



H = -i^^/gmie 



(34) 



The condition ( 12 ) that and cf) represent a physical 
solution to the Bogolubov equation ^ translates into 
reality conditions for hi and 9i. It is also interesting to 
note that the commutation relation for these two fields 



IS 



[ni(x, t), 6'i(x', i)] — ighd{x 



(35) 



as a direct consequence of the commutation relation ([8 
for (j). Thus, hi and 9i are canonically conjugate fields. 

A. Generalized Klein-Gordon equation 



We can now combine the two equations (30 1 and (31 ) 



for 9i and hi in order to obtain a second order (in time) 
differential equation for 9i. More explicitly, we can ob- 
tain hi as a function of 9i from Eq. (31) by formally 
inverting the operator (1 — 8): 



ni 



-Widt + V • W)9i 



(36) 



where 



As an alternative to the approach followed in the pre- 
vious section, we can linearize the continuity ^ and 



w = (i-e)-i = 1]®" 



(37) 
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Note that generally this is a well-defined procedure be- 
cause of the negativity of the operator Q. Indeed, for ho- 
mogeneous profiles, O is proportional to the Laplacian, 
which is obviously a negative operator. As long as the 
profile is sufficiently smooth such that c is slowly varying 
on the healing length scale, this negative character will 
not be altered. Relevant departures from this behaviour 
would require profiles whose density varies significantly 
within length scales comparable with the healing length, 
which is not only easily avoidable in practice but might 
actually be very hard to realize. Keeping these comments 
in mind, we can insert Eq. (36 1 into Eq. (30) to obtain 
the single equation 



This means that an effective curved spacetime geome- 
try is recovered only for condensate perturbations such 
that W is very close to 1. Furthermore, as we will see 
below, the dispersion relation for this W-Klein-Gordon 
equation is the same as that for the Klein-Gordon equa- 
tion modified with fourth-order spatial derivatives. Here, 
however, the Klein-Gordon inner product is modified in 
each t = constant slice by the action of the operator W. 
The Klein-Gordon equation modified with fourth-order 
spatial derivatives, on the other hand, shares the Klein- 



Gordon product in each t =constant slice (41 1 with the 
proper Klein-Gordon equation [16 . 

In spite of these modifications, the W-Klein- Gordon 
product ( 42 ) shares the following properties with the 



[dt + V ■ (v*)]W(9t + V • V)6ii + V • (c^V6'i) .(38) standard Klein-Gordon product: 



This equation is a higher order differential equation 
which generalizes the Klein-Gordon equation. 

In the limit W 1, attained when the gradients in 
9i are relevant only for length scales much larger than 
the healing length, a proper Klein-Gordon equation is 
recovered. Indeed, in this case, this evolution equation 
can be written in the form W 



where is the acoustic metric 



-(c2-v2) 

V 



1 



(39) 



(40) 



The corresponding Klein-Gordon inner product can be 
written as 

= i J d^x^?i*f9t + v- vj^'i , (41) 

where = c^/(^~i) is the determinant of the metric in 
the spatial slice i=constant, n^^ — c~^/(-°~^)(l, v) is its 
normal, and i9„ — n^d^. 



• It is conserved in the lab time t. However, unlike 
the Klein-Gordon product, which is conserved in 
any inertial time, the W-Klein-Gordon product is 
conserved only in the lab time. This is a logical 
consequence of the fact that local Lorentz invari- 
ance is no longer in operation. 



• It is hermitian, i.e. 



^l|^'l)vV-KG 



'llfl/W-KG- 



• It is antilinear in the first argument and linear in 
the second. 



It is not positive definite. Indeed, 

(^ll^i*)W-KG = -(^''ll^l)w-KG 



(43) 



• There exists at least one basis {On} of orthonormal 
solutions such that 



(44) 

The norm of any solution 6i = '^^{aiOu + b*6l^) is 
(0i|ei)w-KG = ^(|a.p-|6.n . (45) 



{0li\9lj}w-KG 


— ^ij , 


{Oli\Olj)w-KG 


= , 


(^ljl^'lj>W-KG 


= 



B. Generalized Klein- Gordon product 

In the general case, we can introduce an inner product 
that generalizes the Klein-Gordon product by taking into 
account that the time derivative term is now modified by 
the operator W. In D spatial dimensions, this W-Klein- 
Gordon inner product turns out to be 



w-KG = ^ / d^x 91 [W{dt + V • V)]9[ . (42) 



It is clear from this expression that it reduces to the 
standard relativistic Klein-Gordon product in the limit 
1. 

It should be stressed that the operator W breaks the 
local Lorentz invariance of the Klein-Gordon equation. 



This norm vanishes for real scalar fields such as the 
physical phase perturbation. 

At this stage, it is worth noting that the W-Klein- 
Gordon inner product has the same properties that we 
stated above for the Bogolubov inner product. In fact we 
will see that these two products are indeed equivalent. 



C. Mode expansion 

We will now discuss the form of the modes and their 
normalization. For the sake of simplicity and definite- 
ness, let us concentrate on a case of particular interest, 
namely, the case of a background profile which becomes 
time-independent in the asymptotic future. Then we 
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can try an ansatz of the form ^e~*""^*'''\ which in the 
pure one-dimensional Klein-Gordon case leads to exact 
orthonormal modes (2c|A:|)~^/^e~*('^*~'^^^. 
Let us define 

<D(t,x) = Ludtu , k(<,x) = —LdWu . (46) 

The condition that the profile becomes stationary in the 
asymptotic future implies that ui{t —i' oo,x.) — ui and that 
k(t,x) is time-independent in this limit. Then Ae"*"" 
is an approximate solution which can be found in the 
regime where k, a), and A are slowly varying functions 
(in space and time). Introduction of this approximate 
solution into the modified Klein-Gordon equation yields, 
to lowest order, the dispersion relation 

(lj - V • k)2 = c^k^rl , (47) 



the Bogolubov equation and use this expression in the 
complex conjugate of this same equation, we obtain 

{[h{dt+^-V}~ iTp] 1 [h{dt + V • V) + iTp] 

+ 2mTp^ip^0, (52) 

where (p = (jj/ipo and Tp = — V • (c^V*). It might 
seem at first sight that the number of degrees of freedom 
is equivalent to that of a complex scalar field (i.e., two 
sets of complex Fourier coefficients). However ip and ip* 
are not independent but are related by the Bogolubov 
equation itself, leaving — as before — ^just one field degree 
of freedom (a single set of complex Fourier coefhcients) 
and, in fact, the scalar product in the space of solutions 
of this equation is the Bogolubov one. 



and hence the form of u{t, x), where 

rk = +v/rTFkV4 



(48) 



The next-to-lowest order yields the prefactor A so that a 
complete set of approximate modes is given by 



7lk = 



(49) 



The next order provides information about the spread 
of these modes in the directions perpendicular to that 
of group propagation. The relativistic limit W ^ 1 ob- 
viously corresponds to Fk 1, i.e., to the long wave- 
length limit ^|k| <^ 1, which gives rise to the well-known 
Klein-Gordon modes. These approximate modes are or- 
thonormal (in the same level of approximation) in the 
W-Klein-Gordon product ( 42 ) as can be seen by evalu- 
ating it at < — > 00. Evaluated at any other finite time 
t, small deviations which are consistent with the level of 
approximation that we are using may appear. 

As mentioned above, any real solution 0i can be writ- 
ten as a linear combination of these modes 



9i — y~^(ak6*ik 



*k'^lk7 



(50) 



and therefore has zero norm. 

For completeness, it is also straightforward to check 
from Eq. ([36| that 



.c\k\ 
Tk 



'Ik 



(51) 



Finally, before examining the relation between the Bo- 
golubov and the Klein-Gordon formalisms, we mention 
the existence of a third alternative approach |17j , which 
leads to a genuinely complex differential equation of sec- 
ond order in time and fourth order in space for the field 
perturbation (f>. Indeed, if we write (f)* in terms of 4> from 



IV. EQUIVALENCE BETWEEN BOTH 
PRODUCTS 



The following question of interest regards the rela- 
tion between the Bogolubov product and the W-Klein- 
Gordon product introduced in the previous sections. 
More specifically, what is the relation between the con- 
cepts of quasiparticle and antiquasiparticle in both for- 
malisms? 

Let us consider two spinors <I> and <&', solutions to the 
Bogolubov evolution equation (10). Then, taking into 



account the relation (33) between both representations, 
i.e., between the perturbation fields (f), 4> and hi, 9i, we 
obtain 

($1$% = \j d^x(0>'-0*0') 



2gh 



(53) 



Finally, the relation (36) between the density and 



phase perturbation allows us to write this Bogolubov 
product as 



2gh 
1 

2^ 



d^:E( 



W{dt 



V)] 



W-KG • 



(54) 
(55) 



So, we see that the Bogolubov and the W-Klein- Gordon 
products are indeed equivalent. 

Also, provided a set of orthonormal modes 6'ik of the 
W-Klein-Gordon equation, it is straightforward to con- 
struct an orthonormal set of modes for the Bogolubov 
equation (10) by means of the relations (33) and (51): 



Wk = e^''"/" 



Wk 



-ieo/h 




(56) 
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Therefore positive (resp. negative) norm modes in the 
W-Klein-Gordon product are mapped to positive (resp. 
negative) norm modes in the Bogolubov product and vice 
versa. This means that, whether we choose to analyze 
quasiparticlc creation processes (e.g., Hawking radiation 
in a black hole configuration) in a modified relativistic 
framework such as in '18] or a condensed-matter context 
such as in 17J, the results should coincide. The reason 
for this coincidence is that, as we have seen, there is 
a one-to-one relation between both formalisms, and the 
concepts of quasiparticle and antiquasiparticle, as well 
as the ambiguities inherent in these definitions (which 
are related to the observer-dependence of the concepts 
involved), are the same in both formalisms. 



V. ON THE CHOICE OF VACUUM STATE 

One can now proceed with the quantization of the sys- 
tem following the standard Fock procedure. One only 
needs to find an orthonormal mode basis to expand the 
field operators that characterize the quantum perturba- 
tions. In this manner, one can define creation and annihi- 
lation operators and a vacuum state for the system. The 
fact that the inner product is not positive definite tells us 
that the selection of a specific set of positive (negative) 
norm modes can be done in many different ways. To il- 
lustrate the procedure, let us consider a one-dimensional 
stationary fiow in the condensate which simulates the 
presence of a black hole, but with the internal singularity 
substituted by a second asymptotic region |19| . Among 
the different sets of modes that can be selected, there are 
three of special relevance, which we will call the "in" set, 
the "out" set and the "stationary" set, whose meaning 
will be discussed shortly. Each positive energy mode of 
each set can be characterized by its frequency (which is 
invariant due to the stationarity of the system) and by 
a discrete label which, depending on the frequency, can 
acquire the values 1, 2 or 1, 2, 3 (see the discussion in [20] 
regarding the number of normalizable modes in different 
configurations) . 

We will now classify these modes, using the following 
subscripts. The p and / subscripts represent "past" and 
"future" (in a scattering process) while u and w repre- 
sent the "right-going" and "left-going" character of wave 
packets centered around the particular mode frequency 
with respect to the lab. 

Let us consider the dispersion relation ( |47[ ) formally 
as a relation between w, k and x (through v[x)). Then, 
for any given real frequency uj and any point x in the 
configuration, the dispersion relation may have two or 
four real roots fc, see Fig. [l] Actually, there exists a 
critical frequency Wc such that, for w > cJj,, independently 
of X, there always exist two real roots only, one positive 
and one negative. We denote by the positive root 
when X — >■ — oo, by k^p the absolute value of the negative 
root when x — -l-oo, by k^f the positive root when x — 
-|-oo and finally by k^f the absolute value of the negative 




FIG. 1: Dispersion relation \41\ for various values of uj and of 
\v\, scaled at c = 1. The intersection points between the line 
UJ — vk (in blue) and the various branches of the dispersion 
relation (in black) mark the real (normalisable) mode solu- 
tions for a given oj. The critical frequency uuc represents the 
frequency at which additional, 'extraordinary' roots appear 
for the maximal |«| attained in the configuration. Then, for 
oj > cjc, there are always only two normalisable solutions. For 
(jj < ujc, there can be either two or four normalisable solutions, 
depending on the value of 



root when x — > — oo. For < w < cj^, depending on x, 
there can be either two or four roots (there is a critical 
position at which there are just three solutions; here we 
will not discuss this critical situation). Two of the roots 
are always equivalent to the previous ones; we will use for 
them the same notation. In case of the existence of four 
real roots, the two additional ones always correspond to 
negative values of k. We will call these additional modes 
'extraordinary' because they are absent in the subsonic 
regime. Let us then denote by k^p the absolute value 
of the most negative extraordinary root when a; — >■ — oo 
and by k^f the absolute value of the second most negative 
extraordinary root when x — > — oo. The corresponding 
wave packets are right-going (kep) and left-going {kef), 
respectively. 

We summarize this classification in table H] 
Given the equivalence displayed in the previous sec- 
tion, we are free to describe the perturbations either by 
a spinor in the Bogolubov formulation or by the phase 
perturbation operator 6 in the Klein-Gordon formulation. 
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t — )■ — oo 


t — > +00 


— > 




— > 




X — )■ H-co [v < c) 




kwp 


kuf 




X — ^ — oo {v > c) and uj > ljc 








kwf 


X — )■ — oo {v > c) and lj < cjc 


kup 

kep 






kwf 



TABLE I: Classification of normalizable mode solutions (real 
roots) in the asymptotic regions. The arrows indicate the left- 
or right-going character of the corresponding wave packets 
with respect to the lab. 
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^ ^ 

<f • • "^^y^ 









For notational simplicity, we omit the subscript 1 for the 
perturbations in this section. 



A. "In" vacuum state 

It can be shown [20] that the positive energy modes 
that constitute the "in" set (such that, in the "in" vac- 
uum defined by them, there are no quasiparticles in the 
asymptotic past) can be described as follows, up to a 
global mode-dependent normalization constant: 



• For oj > Wc, 



a:— >-+oo 



in ^ik^fX 



LJ,2 ^ ^ \^ 



(57) 



(58) 



FIG. 2: "In" basis, u > ujc'- (1) represents a scattering pro- 
cess in which there is an incoming wave from the left; (2) 
represents a scattering process in which there is an incoming 
wave from the right. < u < ujc'. {!') represents a scattering 
process in which there is an incoming wave from the left; (2') 
represents a scattering process in which there is an incoming 
wave from the right; (3') represents a scattering process in 
which there is an extraordinary incoming wave from the left. 
Solid arrows represent ordinary modes while dashed arrows 
represent extraordinary modes. 



right. For < uj < ujc there are three independent and 
orthonormal scattering processes: (1') an incoming wave 
from the left, (2') an incoming wave from the right, and 
(3') an incoming extraordinary wave from the left. All 
these three processes end up being a combination of three 
outgoing waves at future infinity, one of them extraordi- 
nary. Notice that uj in the exponential of the mode 6*^3 
appears with a positive sign, contrarily to the negative 
sign in the other modes. These signs are chosen in order 
to define positive energy (norm) modes. ^ 



• For io < LUc, 



\i\^ — ik.uj fX 



R" 



-ike 



— ik-u 



— ikp 



T e*"* (e-''== 



(59) 



(60) 



(61) 



where the i?'s and T's are mode-dependent reflection and 
transmission coefficients. In other words, every mode be- 
haves in the specific form displayed by the above formulas 
in the appropriate asymptotic limits. These modes can 
be pictorially described as representing elementary scat- 
tering processes (see Fig. |2]). For u} > Wc there are two 
independent and orthonormal scattering processes: (1) 
represents the scattering of an incoming wave from the 
left and (2) the scattering of an incoming wave from the 



B. "Out" vacuum state 

In an equivalent way, one can construct the positive en- 
ergy "out" basis (such that, in the "out" vacuum defined 
by them, there are no quasiparticles in the asymptotic fu- 
ture) |20| . up to a global mode-dependent normalization 
constant: 

• For Lo > uJc, 



nout x^+po ~iu)t 



out ^—iku. 



QOUt 



nout x—^—po —iijt ( ^ — 



9 out a;— >+oo —-itot ( rpont —ik^^x ^ 
9 — > G 1,. e 



^^-ik^fx 

(rriOUt —ikyjpX \ 
J 



iX, , c- 



In the notation of Ref. 1201 our ^J^g corresponds to </:^^ 



(62) 



(63) 
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UJ < U)c i 
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(1') 
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(2') 
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."^ UJ < UJc i 

(3') 
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.V 

(1') 

a; 






y-^ <j » 










(2) 



(2') 





(3') 




X 


<J 







FIG. 3: "Out" basis. Interpretation as in (time-reversed) 
Fig-i 



• For < w < Wc, 

gout g-^ct ^^^k^fx _^ 

QOUt 



(64) 



QOUt ^ — iujt 



fX _|_ ^OUt ^ikupX _|_ J^OUt^lkep 



Jt ^ 



(65) 



gout c*^* fc'^^sfX 



(66) 



FIG. 4: Preliminary "stationary" basis, cj > tOc'- (1) repre- 
sents a scattering process in which there is an outgoing wave 
to the right; (2) represents a scattering process in which there 
is an incoming wave from the right. < uj < ujc'. {!') repre- 
sents a scattering process in which there is an outgoing wave 
to the right; (2') represents a scattering process in which there 
is an incoming wave from the right; (3') represents a scatter- 
ing process in which there are neither incoming nor outgoing 
waves on the right. These modes are not orthogonal. The 
"stationary" orthonormal basis is obtained by linear combina- 
tions of these modes. Solid arrows represent ordinary modes 
while dashed arrows represent extraordinary modes. 



and the additional mode 9^ ^ given by 



-ike fx\ 



The pictorial representations of these elementary scat- 
tering processes can be seen in Fig. [3] They correspond 
to the elementary scattering processes of Fig. [2] running 
backwards in time (note that for this identification we 
have to exchange labels 1 and 2). 

For the case oj < ojc, it is not possible to express 
the positive norm mode 0°"* only in terms of the pos- 
itive norm modes 6*2^]^, 0™2j ^™3- has to use the 
negative norm modes ^L" i ; ^^.2 7 3 ^ well. Therefore, 
the "in" vacuum state defined through the requirement 
'^™i|Oin) = contains quasiparticles coming out at the 
right asymptotic region. This phenomenon is usually 
called mode mixing and has been thoroughly analyzed 
in m. 



C. "Stationary" vacuum state 

Let us finally describe the "stationary" set. As a pre- 
liminary stage, let us choose 



• For UJ > UJc, 



flSt-p 



/lOUt 



and 



gst-p _ 



1,2 I 



• For < UJ < UJc, 

/)St-p nout 



nst-p nin 



(67) 
(68) 



„st-p an— f-Kx) _ 



(69) 



if it has positive norm; or by 



QSt-p X- 



'iku 



{A^e 



B,. 



JkefX\ 



0: 



st-p X — )-+o o 
uj,3 







(70) 



if it occurs that the conjugate of the expression in 
Eq. ( 69 1 is the one with positive norm (this will de- 



pend on the specific shape of the velocity profile). 
The ratio Ai^/Ai^^g is fixed by the condition that 
the scattering lead to no transmission in the right 
asymptotic region. Then, for a given A^, B^^ and 
Bi^ c become fixed. Finally the value of A^ itself 
can be fixed by requiring that the mode be nor- 
malized. Fig. |4] shows a pictorial representation of 
these modes. 

The previous three modes are clearly independent and 
normalized but they are not orthogonal. However, start- 
ing from them, it is easy to define three new modes which 
are orthonormal. Consider the transformation 



, (71) 





^^1 ^ 










ast 






^a;,2 
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^^,3 ) 





Afn 


M12 


Ml3 \ 
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/jSt-p 





M22 


M23 












M33 / 
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/jSt-p 
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with 

Mil = N-^'^ ■ (72) 

Ml2 = iV-l/2(p23Pl3 -Pt2)(l - \P25?)~^ ; (73) 

Mi3 = N-^'^{p1^pI^-p\^){\ - \p2z?r' ; (74) 

M22 = (l-|p23n-'/' ; (75) 
A^23 = ~P23(l-b23n~'/' ; (76) 

M33 = 1 ; (77) 
where 

iV = (l-b23n"'x (78) 

(1 - b23p - bl2p - bl3p +Pl2P23Pt3 +P12P23P13) ; 

P..=P;. = (C7lOw-Ko. (79) 

It is easy to check that the new modes are indeed or- 
thonormal. Then the corresponding annihilation opera- 
tors transform with the transposed inverse of the matrix 
M: 









/ 


st-p \ 




<2 


= {M-'Y 




st-p 
"<^,2 


V 


<3 ) 




V 


st-p / 
«<^,3 / 



The important point here is that M is an upper trian- 
gular matrix. Therefore, the matrix {M~^y appearing 
in this equation is lower triangular. Also, M^, which al- 
lows us to write the annihilation operators a^' f as linear 
combinations of a^'j, is a lower triangular matrix. This 

implies that a^*~f , a|j^~|' depend linearly on a^* i, 1^*2 only, 
and vice versa, never mixing with the third mode, which 
would lead to particle presence in the right asymptotic 
region. 

Consider the "stationary" vacuum state defined by the 
requirement 

<.|Ost>=0. (81) 

The important point is that this state does not contain 
any quasiparticles coming in or out of the right asymp- 
totic region. To realize that this is indeed so, one has to 
check that the Bogolubov /3 coefficients defined by the 
products {e^ieZ'^l *)w-KG and (^L* JC,2*)w-kg are identi- 
cally zero {i = 1,2,3). Given the relation between the 
and the dl^~[ modes, this follows straightforwardly 
from the following argument. 

By the definition of 9'^~l and 9^^~2 it is obvious that 

KrKl*)w-Ko = and (C7lC'^2*)w-Kc = for z = 
1, 2. In order to show that (0^*'3 |^J^,2*)w-kg = 0, evaluate 
the spatial integral that defines this inner product at i — >■ 
—00. Then, it can be seen from the scattering diagrams 
that the involved modes do not have intersecting support, 
and therefore their product vanishes. A similar argument 
applies for (^^*'3 16*°^! *)w-kg by evaluating it at t = -l-oo. 
(Note that this can be seen even more clearly by using a 
wave packet basis.) 



Therefore in the "stationary" vacuum state there are 
no quasiparticles coming in or out from the right asymp- 
totic region. In other words, this state is a generaliza- 
tion to this dispersive theory of the Boulware station- 
ary state for black holes and static stars. In standard 
general relativity, the Boulware state is not regular at 
the horizon. The dispersive nature of Bogolubov theory, 
however, makes this state perfectly regular even at the 
horizon. Consequently, it is a perfectly attainable state 
in such a dispersive theory. 



D. State preparation 

Given a stationary configuration, one could prepare the 
system to be in any of the vacuum states described above. 
The "in" vacuum contains quasiparticles coming out at 
the right infinity {x — >■ -l-oo), while the stationary vacuum 
state does not (it docs not contain quasiparticles coming 
in from the right infinity cither) . In standard general rel- 
ativity, one can prove that starting from a Minkowskian 
spacetime at past infinity and a relativistic field initially 
in its Minkowski vacuum, if one dynamically produces a 
black hole horizon, then the vacuum state at future infin- 
ity is indistinguishable from the Unruh state ^21t i22 . In 
particular, if one takes a Lorentz invariant theory and en- 
gineers, by using external means, a quasistatic collapse 
towards the formation of a black hole, then before the 
formation of the horizon, one can always make the col- 
lapse slow enough so that an initial Minkowskian state 
first acquires and then maintains at each instant of time 
the Boulware vacuum structure, which is the appropriate 
vacuum state for a stationary star. However, no matter 
how slow the collapse is, the quasistatic approximation 
breaks down when the horizon forms due to the infinite 
slow-down of clock rates at the horizon. As a conse- 
quence, the Boulware state is unavoidably modified and 
becomes, after some transient phase, the Unruh state 
with its associated Hawking emission. Note that this is 
also necessary for consistency: the renormalized stress- 
energy tensor in the Boulware state is divergent at the 
horizon and the state itself is not well defined there, be- 
cause of the infinite blue-shift of the corresponding modes 
caused by the presence of the (non-dispersive) horizon, as 
pointed out above. The Unruh state, on the contrary, is 
perfectly regular at the (future) horizon (we don't men- 
tion here states, such as the Hartle-Hawking state, which 
are not vacuum states in the past). Thus, within stan- 
dard general relativity a field in a black hole spacetime 
can only be in the Unruh vacuum state (obviously, mod- 
ulus the presence of any finite amount of particles) and 
never in the Boulware vacuum state. 

In dispersive theories, the previous arguments are no 
longer valid and, depending on the specific preparation 
of the configuration — for instance the specific dynamical 
way in which one sets up the subsonic to supersonic tran- 
sition in the Bose-Einstein condensate fiow — one could in 
principle end up having different final states for the final 
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stationary flow. Indeed, in a Bose-Einstein condensate 
there is no infinite delay of the clocks and so it should 
be possible to maintain the Boulware vacuum structure 
even after the horizon has been formed. In principle, 
one could for example end up in the stationary vacuum 
state described above, which is a regular extrapolation 
of the Boulware vacuum to the dispersive theory. There- 
fore, although in general the formation of a horizon in a 
Bose-Einstein condensate leads to a Hawking-like radia- 
tion associated to the "in" vacuum (as shown numerically 
in |10) and studied in detail in |17|). it remains to be seen 
which the precise sufficient conditions are in order to re- 
cover the presence of such Hawking radiation in systems 
with superluminal dispersion relations. This could for 
instance be relevant when analysing the effect of high- 
frequency superluminal dispersion in proper black hole 
configurations and the influence of the black hole's inter- 
nal region. 



This issue has recently created some controversy. In 
Ref. [18 the present authors argued that, under the as- 
sumption that a quasi-static condition for the creation of 
a black hole applies, one does not need to take into ac- 
count all the modes of the system to calculate the quasi- 
particle content at right infinity, but only those that can 
be traced as rays escaping from the black hole configu- 
ration just before the actual horizon was formed. Quasi- 
particle production due to the formation of a black hole 
in the lab would then appear as just a transient regime 
and disappears in time: the system settles down to the 
stationary state described above. However, other authors 
claim that the relevant state for the final configuration 
is the "in" vacuum state (see for example ^20^). This 
vacuum state produces an ever-lasting stationary stream 
of quasiparticles travelling towards the right asymptotic 
region, mimicking in this way a Hawking flux (the pre- 
cise spectrum acquires some deviations with respect to a 
perfect black body but, in normal situations, these de- 
viations arc relegated to the high-frequency tail of the 
spectrum). If the "in" vacuum state were the only one 
available for the configuration, one could conclude that it 
is not possible to create a stable sub-to-supersonic tran- 
sition in the lab: one would have to maintain it there 
by external means or it would eventually dissolve due to 
backreaction. This is what happens in standard general 
relativity: the existence of stationary black holes is not 
semiclassically consistent, they have to evaporate. 



In the Bose-Einstein-condensate dispersive theory, 
the availability of the "stationary" vacuum implies 
that, in principle, it should be possible to produce a 
semiclassically-stable analogue of a stationary black hole 
in the lab. Which is the precise vacuum state selected by 
the dynamical formation of the black hole horizon con- 
figuration under discussion remains to be calculated. 



VI. CONCLUSIONS AND COMMENTS 

When establishing a gravitational analogy in con- 
densed matter systems, a low-energy Klein-Gordon dy- 
namics emerges for the perturbations around a back- 
ground configuration. In this paper, we have addressed 
the issue of the appropriate formulation for dealing with 
the quantum dynamics of such fluctuations: the one ob- 
tained by directly dealing with the quantum creation 
and annihilation operators for the excited quasiparticles 
(the Bogolubov formalism), or the relativistic approach 
consisting in obtaining a generalized Klein-Gordon equa- 
tion and dealing with it in the relativistic quantum- 
field-theoretic way. This question applies in particular 
to the quasiparticle Hawking radiation when the back- 
ground configuration supports an acoustic black hole, ei- 
ther stationary or externally generated. As we have seen, 
both methods are entirely equivalent, leading to the same 
quasiparticle concept, and hence to the same description. 

It should be stressed that these results are actually 
more general than presented here. Indeed, if instead of 
a Bose-Einstein condensate we had considered an arbi- 
trary barotropic, inviscid and irrotational fluid described 
by an arbitrary enthalpy function h(n), we would have 
drawn the same conclusions. Furthermore, for regimes 
in which the quantum potential is not signiflcant (long 
wavelengths) , the fluid is described in pure hydrodynam- 
ical terms and a proper Klein-Gordon equation is recov- 
ered. Therefore, we can conclude that the relativistic 
analogy with a hydrodynamic fluid is actually more than 
just an analogy: there is a complete equivalence leading 
to a hydrodynamical description of relativistic massless 
scalar fields and vice versa, even at the quantum level 
and in strong gravity regimes such as black holes. 

Once the inner product is defined one can find a com- 
plete and orthonormal set of positive norm modes, and 
an associated set of negative-norm ones. One can then 
define quasiparticle creation and annihilation operators 
and make a choice of vacuum state. This procedure, as 
is well known, can be carried out in many equally valid 
ways, related through Bogolubov transformations. To 
illustrate the procedure, we considered a simple config- 
uration of special interest with regard to the possible 
experimental detection of analogue Hawking radiation: 
a one-dimensional configuration simulating the presence 
of a black hole horizon for acoustic perturbations or, in 
other words, a flow with a subsonic-to-supersonic transi- 
tion. 

In the last section, we have shown that there are several 
sets of modes and vacuum states with a particularly sim- 
ple interpretation: the "in" modes and the "out" modes 
(with their associated "in" vacuum and "out" vacuum 
states), already present in previous analyses in the liter- 
ature [50], and the "stationary" modes and "stationary" 
vacuum state, presented here for the first time. Contrar- 
ily to standard general relativity, the dispersive character 
of the Bogolubov theory allows the existence of a regu- 
lar vacuum state which does not contain any incoming 
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nor outgoing quasiparticles in the external asymptotic 
region, namely, the "stationary" state. This state is a 
regular generalization of the notion of Boulware state in 
black hole physics. Since it is regular, this state can in 
principle be attained by adequately preparing the system. 
In general, the formation of a horizon in a Bose-Einstein 
condensate should lead to a Hawking-like radiation as- 
sociated to the "in" vacuum [Tni H?]. However, in the 
light of the present discussion, it should at least theo- 
retically be possible to produce a black hole analogue 
in a Bose-Einstein condensate without causing the emis- 
sion of a stationary Hawking flux. Whether a black hole 
analogue in a Bose-Einstein condensate radiates or not 
might depend on the specific path followed in setting up 
the configuration. This could be an important issue for 
experimentalists trying to reproduce Hawking radiation 
in a lab. From a relativistic point of view, it could also 
be of considerable importance when analysing the effect 



of modified dispersion at high energies in gravitational 
black hole configurations in which the characteristics of 
the internal region are uncertain. We will return to these 
issues in future studies. 
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